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ABSTRACT 


A  method  of  relative  magnitude  determination  using 
all  available  data  for  a  limited  source  area  is  pvt josed. 
The  attenuation  factor  B,  event  terms  Fj  ,  and  station 
terms  for  the  linear  model 


C.  +  F. 
i  3 


+  H 


are  solved  simultaneous  1>  by  an  iterative  leas'"  - 
squares  procedure;  here  is  the  logarithm  /  *.  the 
amplitude  and  is  the  logarithm  of  the  distance.  The 
method  is  analogous  to  joint  epicenter  deternunation, 
and  so  is  termed  "joint  magnitude  determination"  (JMD) . 
Since  observed  amplitudes,  either  body  ware  o '•*  surface 
wave,  seldom  are  seen  to  follow  a  straight-line  decay, 
the  model  is  a  gross  simplification  in  which  the  station 
terms  actually  operate  somewhat  as  distance- correction 
terms.  The  JMD  method  was  applied  to  Seismic  Data 
Laboratory  (SDL)  data  from  North  American  stations  for 
Nevada  Test  Site  explosions.  The  standard  deviation  of 
station  magnitudes  (Ms  and  m^)  computed  by  JMD  for  an 
event  typically  was  one-half  of  the  standard  deviation 
of  magnitudes  computed  by  routine  methods.  This  reduc¬ 
tion  is  wholly  due  to  the  addition  of  station  terms. 
Surprisingly,  no  consequent  reduction  in  scatter  of  the 
explosion  M  vs  m^  plot  was  attained  by  using  JMD  esti¬ 
mates  rather  than  routine  ones,  and  it  is  concluded 
that  source  parameters  rather  than  bias  due  to  varying 
networks  is  causing  the  scatter.  The  scatter  of  magni¬ 
tude  estimates  is  divided  into  three  physically  meaningful 


parts:  source,  path-station,  and  measurement 0  The 
absolute  value  for  each  part  is  derived  from  the 
seismic  measurements  without  using  yield  information. 
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INTRODUCTION 


Estimation  of  event  magnitude,  either  or  m^,  has 
always  been  hampered  by  the  inability  to  provide  a 
uniform  distance-correction  function  which  will  apply 
satisfactorily  to  all  source  areas.  Especially  at 
regional  distances,  application  of  generalized  distance 
corrections  can  introduce  extreme  bias  into  magnitude 
estimates.  Therefore,  seismologists  have  devised  a 
number  of  disparate  magnitude  formulas  and  distance- 
correction  functions  or  tables  (B&th,  1966  and  1969; 
VESIAC,  1964).  For  example,  in  the  Western  United  States 
detailed  analysis  of  the  earth  structure  and  travel 
times  enabled  Evernden  (1967)  to  formulate  appropriate 
m^  formulas;  and  a  study  of  Rayleigh  wave  amplitudes 
from  NTS  events  by  von  Seggern  (1970)  led  to  a  regional 
distance-correction  function  for  the  Western  United 
States.  In  these  and  similar  studies,  the  object  was  to 
devise  a  set  of  formulas  which  would  cover  a  considerable 
distance  range  so  that  recorded  amplitudes  at  given 
distances  could  be  equalized  by  use  of  the  appropriate 
formula. 

There  are  three  disadvantages  to 'the  above  approach: 
1)  the  determined  formulas  are  tied  to  a  certain  area 
and  cannot  be  reliably  applied  elsewhere  when  either  epi¬ 
centers  or  travel  paths  are  changed;  2)  the  amplitude 
analysis  required  to  establish  the  formulas  can  be  lengthy 
and  cumplicated;  and  3)  the  formulas  can  be  awkward  to 
apply  in  practice  when  the  distance  range  has  been 
segmented.  Tabulated  distance-correction  terms  for 


magnitude  have  nearly  the  same  deficiencies. 

Thus,  we  believe  a  simpler  approach  is  demanded. 

What  we  propose  is  to  establish,  using  the  NTS  data, 

single  amplitude-distance  functions  for  Mg  and  for  m^. 

Rather  than  attempting  to  fit  the  observed  amplitudes 

as  closely  as  possible  with  a  segmented  formula,,  we  will 

invoke  station  terms  to  reduce  the  deviations  from  the 

single  amplitude-distance  functions.  These  station  terms 

should  not  be  considered  to  be  intimately  connected  with 

the  immediate  structure  at  the  station,  but  to  be  the 

cumulative  effect  of  the  entire  travel  path  on  the 

recorded  amplitude  at  a  given  site.  Since  we  employ 

station  terms,  it  may  then  be  argued  that  this  approach 

also  is  limited  to  a  small  source  area  and  certain  propa 

gstion  paths.  However,  this  simple  approach  can  be 

applied  to  area  after  area  with  ease,  new  station  corpec 

tions  and  a  new  amplitude-distance  function  being 

generated  each  time  from  the  available  data.  Since  it 

has  been  established  that  the  M  vs  m,  discriminant 

s  b 

should  be  regionalized  (Liebermann  and  Pomeroy,  1969), 
this  magnitude  determination  approach  is  especially 
convenient  in  monitoring  particular  source  areas  where 
relative  values  of  or  m.  between  areas  need  not  be 
established.  We  are  interested  in  determining  whether 
this  approach  can  produce  Mg  vs  m^  data  which  is  as 
good  as  or  better  than  data  provided  by  use  of  more 
conventional  magnitude  calculations.  The  addition  of 
station  terms  should  reduce  the  variance  of  individual 
magnitude  values  for  an  event  and  should  improve  the 
scatter  in  Mg  vs  m^  plots.  In  this  report,  we  reprocess 
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the  data  used  by  Lambert  (1971)  in  establishing  the  Mg 
vs  relation  for  the  Nevada  Test  Site  in  order  to 
compare  this  new  approach  to  routine  magnitude  estimation. 
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METHOD  OF  JOINT  MAGNITUDE  DETERMINATION 


This  method  is  intended  to  compare  magnitudes  of 
individual  events  in  a  particular  source  area  recorded 
by  several  stations.  The  station  set  need  not  be  the 
same  for  each  event,  but  some  stations  should  be  common 
to  several  of  the  events. 

The  following  description  of  the  method  of  obtain¬ 
ing  the  magnitude  formula  pertains  to  both  Mg  and  m^.  As 
the  amplitude-distance  model  for  both,  we  use  an  expo¬ 
nential  decay  function  which  approximates  the  real  data 
reasonably,  as  will  be  seen  later: 


c. 

1 


f . 
1 


h 


where 


x^.  =  the  calculated  amplitude  at  station  i  for 
event  j , 

c^  =  the  station  factor, 
fj  =  the  event  factor, 

h  =  a  constant  for  a  particular  data  set, 
r^  =  the  distance  to  station  i  from  event  j, 

B  =  a  constant  for  a  particular  data  net. 

Linearize  the  equation  by  taking  the  logarithm  of  both 
sides : 


Rewrite  the  equation  ir  this  simple  manner: 


A  A  A  /\  A 

X.  .  =  C.  +  F..  -  II  -  BR.  . 
il  1  1  il 


where  the  substitutions  are  evident.  Here  . nd  henceforth, 
the  carat  symbol  denotes  an  estimate.  Now  let  Z^.  be  the 
log  of  che  ebserved  amplitude;  then  form  the  error  equation: 


-  Z.  .  -  C. 

ij  i 


II  +  BR.. 


One  could,  at  this  point,  construct  the  normal  equations 
for  this  overdetermined  linear  system  and  solve  explicitly 

A  A  A  /\ 

fox  the  B,  H,  C.’s  and  F.'s  which  minimize  the  sum  of  the 
i  1 

squares  of  (1),  However,  it  was  judged  that  the  matrices 
to  be  inverted  would  be  too  large  (nearly  200  x  200  for 
the  data  used  in  this  report)  to  handle  accurately  or 
efficiently;  and  we  used  the  following  iterative  scheme, 
based  on  analytic  solutions  of  the  least-squares  problem. 
This  scheme  allows  us  to  solve  for  up  to  a  combined  total 

/A  ✓N 

of  200  station  and  event  terms,  plus  B  and  II,  in  less 
than  two  minutes  of  computer  time.  We  desire  to  minimize: 


iPh  =  nw,,cz',.,  -  c4 


ij  ij 


A  a  A  A 

F.  -  II  +  BR. 

1  i  y 


where  W . .  is  a  weighting  factor,  equal  to  one  if  station 
i  recorded  event  j  or  to  zero  if  it  did  not.  By  expanding 
the  right-hand  side  of  (2) ,  taking  derivatives  with  respect 
to  B,  and  se  ting  the  result  equal  to  zero,  we  obtain: 

.  -  wwuhj 

H  22  - « - — - - - - - - um*m - -  J  . 
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Similarly,  we  obtain  for  H: 


By  substituting  equation  (4)  in  (3)!,  we  obtain  an  expres¬ 
sion  for  B  which  depends  on  only  the  given  data  plus 
station  and  event  terms: 


-  cnw±j3  CZXW  Z  R  -  nfJW.1R.i  -  R.  )) 

ii  13  ii  1J  1J  1]  ii  ]i  13  ii  1i  1J  1J 

/N  /S 

To  derive  least-squares  expressions  for  the  F ^  * s  and  C^’s 
we  must  sum  equation  (2)  over  j  or  i  only,  respectively. 
Expanding  '2)  with  the  j  summation,  taking  the  deriva- 

/N 

tives  with  respect  to  F . ,  and  equating  to  zero  gives: 

J 

„  ®|Wi j  Ri j  *  pVij  -  '  |Wijei 

Fj  - pT - - - 

l  13 

/N 

Similarly,  we  obtain  for  C^: 


B?Wi3Rij  *  "  IW, 


Ci  ■ 


13  13 

"W 


13  3 


/\  A  A 


The  B,  H,  Fj*s  and  C^’s  are  not  independent  and 
must  be  arrived  at  iteratively.  One  would  like  to  start 

A 

with  some  initial  approximate  values  for  the  F.’s  and 

J 


mri,^i.T'-^ 


accuracy  of  these  initial  F^’s;  they 


s\  S'* 

C^’s.  However,  in  most  cases  the  C^’s  will  be  unknown  or 
only  crudely  known,  and  in  practice  one  can  safely  start 

A  A 

by  setting  all  C^'s  equal  to  zero.  Initial  F_.'s  may  be 
available  from  routine  magnitude  calculations  at  one  or 
more  stations  for  all  the  events.  In  practice,  it  was 
found  that  the  program  converged  to  equivalent  solutions 
regardless  of  ths 

in  fact  can  also  be  safely  set  to  zero  at  the  start  of 
the  iterative  solution.  We  found  that  the  final  values 

S\  A 

of  the  F.’s  and  of  H  reflected  the  initial  values;  that 
3  /N 

is,  if  all  F.’s  are  initially  overestimated,  the  final 
A  3 

F.'s  will  be  overestimated  with  a  consequent  adjustment 
3  /"N 
in  H. 

/\ 

Equations  (5)  and  (6)  show  that  the  intercept  H  and 

✓\  A 

the  slope  B  are  dependent  on  the  distance-sensitive  C. 

S\  S\  1 

terms.  It  was  found  that  estimation  of  H  and  B  must  be 

/N 

isolated  from  calculation  of  the  C^'s  to  obtain  convergence 
of  the  solutions.  This  is  done  by  using: 

[cnw^R.jciR.z,,  -  - 


B  =  jj  *3  13-^  13  13 


13  3 


Cllw  Jcnw.-Rf.)  -  CHW  R  ) 
ji  13  ji  13  13  j.i  13  13 


(8) 


cnwi;j)aiw  Z  .R  -  &  l\i  R  )] 
ji  13  ji  13  13  13  j  3i  13  13 


„  BHWijRij  *  m-jh) 

H  =  - li - 


ji  J  J 


IN 


(9) 


mjOHm 


The  program  iterates  through  equations  (8) ,  (9) ,  (6) ,  and 
(7),  in  that  order,  until  changes  in  the  solutions  are 

negligible.  A  final  iteration  then  returns  to  equations 

^  /\ 

(5)  and  (4)  to  estimate  B  and  H;  the  value  of  the  addi-’ 
tional  terms  in  (8)  and  (9)  as  compared  to  (5)  and  (4) 

A  A 

are  then  relatively  small  and  do  not  really  change  B  or  H 
significantly. 

We  now  point  out  the  analogy  between  this  method  for 
obtaining  station  magnitude  terms  together  with  source 
magnitudes  and  the  Joint  Epicenter  Determination  method 
of  Douglas  (1967)  for  simultaneously  determining  station 
corrections  to  travel  time  and  epicenter  corrections  to 
preliminary  estimates.  Because  of  the  similarity,  we  refer 
to  our  method  as  Joint  Magnitude  Determination  (JMD) . 
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APPLICATION  TO  NTS  COMPRESS I ONAL-WAVE 
AND  RAYLEIGH-WAVE  DATA 

The  staff  of  SDL  has  analyzed  and  tabulated  seismic 
data  on  79  NTS  underground  shots  (FAULTLESS  excluded) 
recorded  at  over  300  stations.  We  will  employ  the  JMD 
program  on  the  compressional-wave  (P„  and  P)  and 
Rayleigh-wave  amplitudes  and  solve  for  B,  H,  F.’s  and 

/V  J 

C^’s  in  both  cases.  Note  that  by  "amplitude"  we  mean 
the  particle  displacement  (A)  and  not  the  particle 
velocity  (A/T)  which  is  often  referred  to  as  "amplitude" 
and  routinely  used  in  magnitude  calculations  at  the  SDL 
and  elsewhere.  It  was  shown  by  von  Seggern  (1970)  that 
the  displacement  measure  had  a  smaller  variance  than 
the  velocity  measure  at  least  for  Rayleigh  waves.  We 
will  show  this  again  with  results  of  JMD;  and  also  in 
this  section  we  will  investigate  the  possibility  of 
using  the  maximum  trace  amplitude  divided  by  the  system 
magnification  at  the  calibration  period  (A),  regardless 
of  observed  period,  another  measure  considered  by  von 
Seggern  (1970). 

Rayleigh-wave  data 

Only  45  of  the  79  shots  had  reported  Rayleigh-wave 
amplitudes.  There  were  712  measurements  available  to 
begin  with.  Through  previous  experience  with  the  data, 
we  felt  that  many  of  the  measurements,  while  not 
anomalous  in  a  real  sense,  are  in  fact  erroneous  due  to 
calibration  errors,  analysts'  error,  transcription  or 
keypunch  errors,  etc.  It  was  therefore  desirable  to 


-9- 


aSSB&asfttttoai^^ 


eliminate  the  tails  of  the  data  distribution.  First.  M 

7  s 

was  computed  in  the  routine  manner  using  our  formula 
(von  Seggern,  1970)  for  A<15°  and  Gutenberg’s  formula 
for  A>15°.  The  average  magnitude  for  each  event  was 
determined,  and  then  all  measurements  which  produced  a 
station  magnitude  more  than  one  unit  from  the  average 
were  deleted.  Fifteen  measurements  were  discarded  in 
this  way.  Three  additional  measurements  were  discarded 
because  they  essentially  duplicated  another  site  for  an 
event  (e.g.,  ST1TX  discarded  while  ST2TX  retained).  A 
different  situation  occurred  with  LRSM  stations  which 
were  repositioned  during  their  recording  career;  such 
stations  are  considered  to  be  the  same  (e.g.,  HL-ID 
and  HL-2ID) ;  and  all  the  measurements  were  retained  under 
one  station  name.  Two  more  measurements  were  discarded 
when  two  events  (both  having  just  one  reported  amplitude) 
were  re-examined:  (1)  DES  MOINES  -  an  excursion  at  KN-UT 
which  was  apparently  picked  originally  as  a  Rayleigh  wave 
is  spurious  because  its  period  and  travel  time  do  not 
coincide  with  other  shots  and  (2)  MINK  -  a  Rayleigh  wave 
picked  at  HL-ID  must  really  be  from  the  Vancouver  earth¬ 
quake  of  about  the  same  origin  time.  After  all  these 
deletions,  692  measurements  remained  for  43  events.  Pro¬ 
gram  JMD  was  then  run  on  this  data.  Initial  input  esti- 

mates  of  the  F.'s  were  taken  from  SDL’s  previous  magnitude 
«/ 

estimates,  except  a  constant  was  subtracted  so  that  they 
would  be  distributed  with  approximately  zero  mean;  and 
the  initial  estimates  of  the  C^'s  were  all  zero.  It  was 
felt  that  some  additional  truncation  of  the  tails  of 
the  data  based  on  their  statistical  distribution  as 
determined  by  the  JMD  program  was  desirable.  The 


-10- 


deviations  of  individual  magnitude  estimates  using 


H 


+ 


(10) 


from  the  term  for  the  j’th  event  appeared  to  be 
normally  distributed;  on  this  assumption  the  standard 
deviation  was  calculated  to  be  0.16.  All  measurements 

A  A 

which  resulted  in  a  value  of  (F.  .  -  -F.)  greater  than 

*  J  J  ^ 

approximately  twice  the  standard  deviation  of  the  F^’s 
(+_  0.35  was  actually  chosen)  were  discarded;  these 
totaled  19  measurements.  At  this  point  six  stations 
outside  North  America,  each  recording  Rayleigh  waves 
from  only  one  event,  were  discarded  from  the  data  base. 
These  were  00-NW,  TE-GL,  FH-PN,  AD-IS,  HW-IS,  and  GG-GR. 

It  is  important  to  mention  that  the  two  island  stations 

/N 

AD-IS  and  HW-IS  had  magnitudes  F^.  (before  station  correc¬ 
tions  were  applied)  of  approximately  0.6  and  0.8  below 

/\ 

their  respective  event  terms  F^ ;  this  anomalous  charac¬ 
ter  may  be  due  to  the  primarily  oceanic  paths  from  NTS. 
After  these  deletions,  667  measurements  were  left  for 
43  events.  Program  JMD  was  run  again,  and  this  time 
three  measurements  fell  over  the  +_  0.35  criterion,  due 
to  slight  readjustments  in  station  and  event  terms. 

These  three  were  discarded,  and  the  final  run  of  JMD' 
produced  no  measurements  over  the  ^  0.35  criterion.  In 
all,  we  have  thus  discarded  approximately  5  per  cent  of 
the  data  because  they  were  anomalous  and  possibly 

/A  A 

erroneous.  The  final  values  of  B  and  H  were  0.90  and 
4.55  respectively,  with  a  95  per  cent  confidence 
interval  of  +  0.03  on  the  slope  B.  The  final  event 
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terms  (equivalent  to  Mg)  are  given  in  Table  I,  and  the 
final  .station  terms  are  given  in  Table  l'l.  Thus,  for 
Rayleigh  waves  from  NTS  to  North  American  stations,  we 
have  this  unified  distance  relation  by  substitution  of 
the  calculated  slope  and  intercept  into  the  original  model 

log  A. .  =  4.55  -  0,90  log  r. .  +  F.  +  C. 

where  A  is  amplitude  in  millimicrons  (peak-to-peak)  and 
r  is  epieentr.al  distance  in  degrees.  In  Figure  la  we  show 
a  plot  of  the  logarithm  of  observed  amplitudes  norma- 

A 

lized  by,  the  event  term,  log(A..)  -  F.,  versus  log  r.. 

J  J  J 

(not  all  664  points  are  shown  since  some  symbols  repre¬ 
sent  more  than  one  data  point  in  this  computer  printout 
plot  routine).  In  Figure  lb  we  show  the  same  data 

A  A 

corrected  by  the  station  term,  log(A.  .)  -  F.  -  C.,  versus 

1 J  J  1 

log  r^y  It  is  evident  that  the  addition  of  the  station 
terms  greatly  decreases  the  scatter  of  the  data. 

The  0.90  exponential  decay  of  the  Rayleigh-wave 
amplitude  data  deserves  discussion,  Gutenberg  (1945) 
suggested  1,66  for  teleseismic  data  (A>15°),  but  Figures 
la  and  lb  do  not  show  a  tendency  for  the  North  American 
data  to  approach  this  value,  even  for  the  teleseismic 
range  from  NTS,  Our  present  result  is  less  than  the 
revised  slope  of  1.09  for  regional  distances  using  the 
same  NTS  Ray'leigh-wave  data  (von  Seggern,  1970);  how¬ 
ever,  that  value  was  arrived  at  by  an  entirely  different 

means.  Our  result  does  agree  very  well  with  that  of 

-  90 

Basham  (1971),  however.  Such  a  slow  decay  as  r  *  can 
only  be  explained  by  a  relatively  high  Q  along  the 
travel  path  of  Rayleigh  waves  at  periods  of  about  12-18 
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TABLE  I 


Rayleigh-wave  Event  Terms 


NUMBER  0I: 


EVENT 

F 

STATIONS 

AUK 

-  ,45 

21 

BENIIAM 

1.60 

7 

BILBY 

.84 

28 

BOURBON 

.17 

11 

BOXCAR 

1.71 

14 

BRONZE 

.15 

27 

BUFF 

-  .08 

20 

CHARCOAL 

-  .34 

17 

CHARTREUSE 

.04 

18 

COMMODORE 

.86 

16 

CORDUROY 

.48 

21 

CUP 

-  .07 

24 

DILUTED  WATERS 

-  .96 

14 

DUMONT 

.41 

18 

DURYEA 

-  .10 

16 

*FISIIER  ' 

.  -1.26 

1 

FORE 

-  .22 

28 

GREELEY 

1,49 

19 

HARDIIAT 

-  .63 

31 

HAYMAKER 

-  .41 

16 

HALFBEAK 

1.11 

14 

KLICKITAT 

-  .24 

25 

KNICKERBOCKER 

.51 

11 

*MADISON 

-  .68 

3 

MARSHMALLOW 

-1.09 

5 

MERRIMEC 

-  .96 

7 

MISSISSIPPI 

-  .32 

34 

NASH 

-  .13 

6 

* PAMPAS 

-  .60 

2 

PALANQUIN 

-1.21 

10 

*PAR 

-1.25 

16 

PILEDRIVER 

.23 

18 

PINSTRIPE 

-  .82 

11 

*REDHOT 

-1.41 

1 

REX 

-  .12 

9 

SCOTCH 

.83 

16 

SCROLL 

-1.23 

5 

*SEDAN 

-  .32 

20 

*STUTZ 

-1.40 

3 

TAN 

.24 

20 

TURF 

-  .28 

31 

WAGTAIL 

-  .10 

19 

WISHBONE 

-1.04 

11 

*Not  used  in  M_  vs  m.  fit 
s  b 


TABLE  II 


Rayleigh-Wave  Station  Terms 


STATION 

AE-NC 

AN-MA 

AR-WS 

AT-NV 

AX2AL 

AY-SD 

AZ-TX 

BE-FL 

BF-CL 

BG-ME 

BL-WV 

BMO 

BM-TX 

BR-PA 

BX-UT 

CG-VA 

CPO 

CP-CI. 

CR-NB 

CU-NV 

CV-TN 

DH-NY 

DR-CO 

DU-OK 

DV-CL 

EB-MT 

EF-TX 

EK-NV 

EN-MO 

EP-TX 

EU2AL 

EY-NV 

FK-CO 

FM-UT 

FR-MA 

FS-AZ 

GD-VA 

GE-AZ 

GI-MA 

GN-NM 

GR2TX 

GV  'TX 

HB-OX 

HH-ND 

HK-WY 

HL2ID 

HN-ME 

HR-AZ 

HT-NM 

HV-MA 

HY-MA 

JE-LA 

JP-AT 

JR-AZ 

KC-MO 

KG-AZ 


NUMBER  OF 
EVENTS 


.08 

-.06 

-.02 

-.29 

-.11 

-.34 

-.20 

.02 

.12 

.17 

-.04 

-.11 

.39 

.42 

.12 

.20 

.17 

.10 

.28 

.03 

-.01 

.24 

-.15 

-.28 

-.41 

-.16 

.07 

-.32 

.12 

.15 

.07 

.27 

.33 

-.25 

-.09 

.02 

-.12 

.10 

-.17 

.21 

.34 

.19 

.09 

-.09 

-.24 

-.07 

-.20 

.15 

-.14 

-.18 

-.03 

.44 

-.02 

.01 

.20 

-.14 


STATION 

KM- CL 

KN-UT 

LAO 

LC-NM 

LG-AZ 

LP-TX 

LS-NH 

LV-LA 

MM-TN 

MN-NV 

MO- ID 

MP-AR 

MU-WA 

MV-CL 

ND-CL 

NL2AZ 

NP-NT 

PG-BC 

PI-KY 

PM-WY 

PT-OR 

RG-SD 

RK-ON 

RT-NM 

RY-ND 

SA4TX 

SE-MN 

SF-AZ 

SG-AZ 

SI-BC 

SJ-TX 

SK-TX 

SN-AZ 

SS-TX 

ST1TX 

SV-AZ 

SV3QB 

SW-MA 

TC-NM 

TFO 

TF-CL 

TK-WA 

TN-CL 

TU-PA 

UBO 

VN-UT 

VO-10 

VT-OR 

WF-MN 

W112YK 

WI-NV 

WMO 

WM-AZ 

WN-SD 

WO-AZ 

WW-UT 


NUMBER  OF 
EVENTS _ 


.02 
-.02 
-.30 
.21 
.09 
.08 
.39 
.22 
.19 
.10 
-.17 
,09 
.33 
.09 
.20 
.20 
.01 
-.10 
-.19 
-.08 
-.05 
-.09 
-.20 
-.15 
-.56 
.33 
.03 
.37 
.02 
.14 
.08 
-.23 
.14 
.04 
.34 
-.45 
-.24 
-.10 
.10 
.01 
-.11 
-.20 
-.14 
.42 
-.25 
-.62 
.18 
- ,  43 
.23 
.12 
.04 
-.21 
.06 
-.10 
.12 
-.29 


Figure  la.  Logarithms  of  observed  amplitudes  minus  event  terms 
(Aij  -  Fj)  versus  log-distance  for  Rayleigh  waves  from  aTS 
explosions  to  North  America  stations. 


OS 


lb.  Same  data  as  Figure  la,  but  with  station  terms. 


seconds  and  by  assuming  that  the  maximum  amplitude  on  the 
recording  is  an  Airy  phase  at  teleseismic  distances. 

Linear  amplitude-distance  functions  were  determined 
using  the  final  664  readings  after  the  data  were  con¬ 
verted  to  A  and  A/T  using  the  reported  periods.  The 
results  are  summarized  in  Table  III  for  comparison  with 
the  amplitude-distance  function  using  A.  Since  the 
confidence  interval  is  proportional  to  the  standard 
deviation  of  the  data,  Table  III  shows  that  A  and  A  are 
as  good  as  or  better  measures  than  A/T,  in  a  statistical 
sense,  when  estimating  magnitudes, 

Compressional-wave  data 

Essentially  the  same  procedure  was  followed  with 
the  compressional-wave  data  for  the  79  NTS  events.  Before 
any  processing,  all  stations  outside  of  the  North 
American  continent  were  deleted.  Event  SHREW,  with  one 
probably  spurious  reading,  was  discarded.  In  addition, 
since  there  were  many  more  stations  having  recorded  P 
waves  than  LR  waves,  it  was  decided  to  retain  only 
those  stations  in  this  analysis  which  had  at  least 
three  compressional-wave  readings.  Several  hundred  obser¬ 
vations  were  discarded  in  this  way.  After  all  these 
deletions,  there  remained  1474  readings  from  106  sta¬ 
tions  for  78  events. 

Magnitude  was  then  calculated  for  each  observation 
in  the  routine  manner,  using  Evernden's  (1967)  regional 
correction  formulas  in  conjunction  with  Gutenberg  and 
Richter’s  (1956)  correction  factors  for  teleseismic 
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distances.  In  order  to  eliminate  possibly  erroneous  data, 
as  for  the  Rayleigh  waves,  all  readings  which  resulted 
in  a  station  magnitude  deviating  more  than  one  unit  from 
the  event  average  were  discarded.  A  total  of  75  readings 
were  discarded  in  this  way,  and  four  more  were  discarded 
because  two  stations  were  left  with  only  two  readings 
each.  Program  JMD  was  run  on  these  remaining  data. 
Normality  of  the  magnitude  residual  distribution  was 
again  assumed,  with  a  calculated  0.22  standard  devia- 

✓V  /\ 

tion;  and  all  measurements  giving  a  value  of  (F^.  “  Pj) 

greater  than  approximately  two  standard  deviations  of 
/\ 

the  's  (+_  0.45  was  actually  used)  were  discarded. 

This  comprised  34  measurements,  and  two  more  were 
discarded  when  one  station  had  only  two  readings  left. 

JMD  was  repeated  on  the  remaining  data,  and  only  three 
readings  failed  to  meet  the  +  0.45  criterion.  In  all, 
approximately  8  per  cent  of  the  data  were  discarded 
because  they  were  anomalous  and  possibly  erroneous.  Using 
the  1356  remaining  readings,  the  final  run  of  JMD  resulted 

A  A 

in  values  for  B  and  H  of  0.87  and  2.11  respectively, 
with  +  0.03  95  per  cent  confidence  limits  on  the  slope  B„ 
The  final  even*  terms  (equivalent  to  m^)  are  given 
in  Table  IV  ,  and  the  final  station  terms  in  Table  V . 
Thus,  for  compressional  waves  from  NTS  to  North  American 
stations  we  have  this  unified  distance  relation  by 
substituting  the  above  calculated  values  of  slope  and 
intercept  into  the  original  model: 

log  A..  =  2.11  -  0.87  log  r..  +  F^  +  C. 

J-j  lj  j  1 

where  A  is  in  millimicrons  (zero-to-peak)  and  r  is  in 
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EVENT 

*AARDVARK 
*ACUSIII 
* AGOUTI 
♦ALLEGHENY 
* ANTLER 

* armadillo 

AUK 

BENIIAM 

BILBY 

*BOBAC 

BOURBON 

BOXCAR 

BRONZE 

BUFF 

♦CASSELMAN 
CHARCOAL 
CHARTREUSE 
♦G11ENA 
* CHINCHILLA 
*CHINCIfiLLA2 
♦CIMARRON 
♦CLEARWATER 
♦CODSAW 
COMMODORE 
CORDUROY 
CUP 

♦DESMOINHS 
DILUTED  WATERS 
♦DOORMOUSE 
♦DOORMOUSE  PRIME 
DUMONT 
DURYEA 
* FEATHER 
♦FISHER 
FORE 
GREELEY 
HARDHAT 

haymaker 

HALFBEAK 


*  Not  used  in  M 


TABLE  IV 

Comp res sional- Wave  Event  Terms 


F 

NUMBER  01 
STATIONS 

.43 

23 

-  .63 

11 

-  .69 

17 

-1.00 

9 

-  .28 

5 

-  .45 

23 

.15 

28 

1.73 

7 

1.08 

30 

-  .67 

14 

.31 

14 

1.72 

12 

.48 

29 

.52 

20 

-  .58 

14 

.33 

19 

.73 

18 

-1.21 

7 

-1.05 

17 

-1.42 

10 

-  .20 

23 

.42 

17 

-  .94 

17 

1.13 

16 

1.06 

21 

.47 

29 

-  .57 

9 

-  .29 

20 

-  .61 

22 

-  .49 

21 

.94 

21 

.49 

20 

-1.79 

6 

-  ,71 

16 

.48 

34 

1.75 

15 

-  .02 

23 

.15 

26 

1.51 

16 

EVENT 

*HYRAX 

*KAWEAH 

klicketat 

KNICKERBOCKER 

*MAI) 

♦MADISON 

marshmallow 

merrimec 

♦MINK 

MISSISSIPPI 

NmSII 

♦packrat 

PALANQUIN 
* PAMPAS 
♦PAR 

* PASSAIC 
*PEBA 

PILEDRIVER 

PINSTRIPE 

♦REDI10T 

REX 

* RINGTAIL 
* ROANOKE 
♦SACRAMENTO 
* SANTEE 
SCOTCH 
SCROLL 
* SEDAN 
■'STILLWATER 
♦STOAT 
* STONES 
♦STUTZ 
TAN 
TURF 
WAGTAIL 
* WICHITA 
WISHBONE 
*YORK 
♦YUBA 


NUMBER  OI: 
E  STATIONS 


-  .64 

15 

-1.00 

14 

.26 

26 

.  66 

13 

-1,77 

5 

-  .44 

16 

-  .28 

16 

-  .52 

14 

-1.30 

6 

.52 

27 

.52 

15 

-  .59 

16 

-  .66 

17 

-  .59 

18 

.01 

26 

-  .53 

21 

-  .48 

19 

.98 

20 

-  .15 

11 

-1.01 

9 

.16 

19 

-1.09 

12 

-1.12 

9 

-  .80 

8 

-  .78 

17 

1.11 

16 

-  .48 

12 

*  .17 

13 

•  .65 

17 

■  .91 

16 

.27 

20 

•  .29 

17 

.97 

21 

.19 

30 

.62 

25 

.74 

14 

.30 

19 

.45 

22 

.61 

21 

vs  mb  fit. 


degrees.  Figures  2a  and  2b  show  for  corapressional  waves, 
as  Figures  In  and  lb  did  for  Rayleigh  waves,  that  the 
addition  of  station  terms  reduces  the  scatter.  Figure  2a 
is  typical  of  the  variation  of  body-wave  amplitude  versus 
distance  plots,  in  that  amplitude  is  not  a  smooth  or 
simple  function  of  distance.  This  character  is  due  to  the 
fact  that  various  branches  of  the  travel-time  curve 
contribute  first  arrivals  out  to  around  20°  distance 
and  that  ray  paths  to  various  distances  sample  widely 
changing  Q  properties  of  the  mantle.  Certainly  the 
straight-line  approximation  calculated  here  is  not  a 
good  fit  to  the  data;  but,  as  Figure  2b  shows,  when 
this  approximation  is  used  in  conjunction  with  station 
terms  which  are  operating  somewhat  as  distance**correction 
terms,  a  small  scatter  in  station  magnitudes  for  any 
given  event  should  result. 

The  statistical  stability  of  the  two  other  measures, 

A  and  A/T,  were  assessed  after  converting  A  to  these 
measures  using  the  reported  periods.  The  JMD  program 
was  run  with  the  final  1356  observations  used  for  the 
A  versus  distance  fit.  Table  III  compares  these  results, 
and  the  confidence  limits  on  the  slope  show  that  all 
three  of  these  measures  for  compressional  waves  have 
approximately  the  same  statistical  stability.  The  A/T 
measure,  however,  does  show  somewhat  less  variance  for 
m^  than  A,  whereas  the  reverse  is  true  for  Mg*  Gutenberg 
(1945)  stated  that  his  body-wave  amplitudes  were  more 
consistent  when  divided  by  the  observed  period,  and  the 
A/T  measure  became  routine  for  m^  calculations  thereafter. 
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Figure  2a.  Logarithms  of  observed  amplitudes  minus  event  terms 
(Aij  -  Fj)  versus  log-distance  for  corapressional  waves  from  NTS 
explosions  to  North  America  stations. 


Distribution  of  magnitude  residuals  and  correlation 
of  station  terms 


Recall  that  we  assumed  a  normal  distribution  for 
both  the  compressional-wave  and  Rayleigh-wave  magnitude 

/N  /N 

residuals  (F^  -  F^).  We  plot  in  Figure  3  the  histograms 
of  the  residuals,  using  the  final  664  and  1356  readings. 
Both  distributions  appear  to  be  normal.  Since  the  data 
were  truncated  beyond  approximately  two  standard  devia¬ 
tions  from  the  mean,  a  slight  modification  of  the 
Smirnov-Kolmogorov  test  for  normality  was  made  before 
it  was  applied  to  this  data.  The  compressional-wave 
residuals  passed  the  normality  test  at  well  above  a 
0.2  level  of  acceptance,  while  the  surface-wave 
residuals  passed  at  only  a  0.05  level  of  acceptance. 
(Levels  were  taken  from  Lilliefors,  1967.)  The  isolated 
peak  around  the  mean  of  the  surface-wave  residuals  in 
Figure  3  undoubtedly  accounts  for  this  low  level  of 
acceptance.  If  the  magnitude  is  regarded  as  normally 
distributed,  then  amplitudes  must  have  a  log-normal 
distribution*  This  log-normal  behavior  has  been  found 
for  two  short-period  data  sets  entirely  different  than 
the  NTS  one  considered  here.  Freedman  (1967)  found  WWNSS 
observations  to  have  this  distribution,  and  Klappenberger 
(1967)  showed  that  LASA  sub-array  amplitudes  for  indi¬ 
vidual  events  had  this  distribution. 

We  checked  the  correlation  of  Rayleigh-wave  and 
compressional-wave  station  terms  at  those  stations 
which  appear  in  both  Table  II  and  V.  The  correlation 
coefficient  w as  -0.10,  which  indicates  that  the  two 
groups  of  station  terms  are  highly  uncorrelated. 
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TABLE  V 


Compressional-Wave  Station  Torms 


STATION 

C 

NUMBER  OF 
EVENTS 

STATION 

C 

NUMBER  OF 
EVENTS 

AN-MA 

-.12 

3 

KN-UT 

.57 

72 

AP-OK 

.31 

4 

LAO 

-.55 

.13 

AR-WS 

-49 

4 

LC-NM 

-.80 

34 

AT-NV 

.75 

16 

LG-AZ 

-.10 

10 

AX2AL 

.33 

10 

LP-TX 

-.14 

5 

Bli-PL 

.19 

7 

MC-SI) 

.04 

3 

BF-CL 

.64 

12 

ML-NM 

-.65 

6 

BL-WV 

-.15 

13 

MM-TN 

-.04 

7 

BMO 

-.79 

30 

MN-NV 

.62 

75 

BM-TX 

-.03 

3 

MO- ID 

-.36 

5 

BP-CL 

.64 

3 

MP-AR 

-.01 

14 

BR-PA 

.04 

8 

MV- CL 

-.37 

34 

BX-UT 

-.13 

3 

ND-CL 

-.14 

4 

CG-VA 

.12 

3 

NG-WS 

.04 

6 

CK-BC 

-.11 

4 

NL2AZ 

-.09 

10 

CN-WS 

.36 

4 

NP-NT 

» 34 

27 

CPO 

.06 

31 

PG-BC 

-.06 

16 

CP-CL 

.12 

41 

PI-WY 

-.17 

5 

CR-NB 

23 

9 

PM-WY 

-.52 

16 

CT-OK 

-.06 

4 

PT-OR 

.02 

25 

CU-NV 

.73 

4  . 

RG-SD 

-.38 

8 

DII-NY 

-.28 

10 

RK-ON 

.70 

28 

DR-CO 

-.44 

37 

RT-NM 

-1.23 

7 

DU-OK 

.40 

4 

RY-ND 

.63 

4 

DV-CL 

.77 

5 

SE-MN 

.19 

17 

EB-MT 

-.02 

5 

SF-AZ 

.07 

13 

EF-TX 

-.48 

6 

SG-AZ 

.16 

8 

EK-NV 

.52 

5 

SI-BC 

-  n  18 

10 

EN-MO 

-.04 

3 

SJ-TX 

.40 

6 

EU-AL 

-.06 

3 

SK-TX 

.02 

4 

EY-NV 

.24 

5 

SN-AZ 

-.13 

10 

FK-CO 

-.20 

6 

SS-TX 

-.38 

15 

FM-UT 

-.32 

36 

ST-NV 

.24 

3 

FO-TX 

-.57 

4 

SV-AZ 

-.28 

4 

FR-MA 

-.02 

5 

SV3QB 

.21 

23 

FS-AZ 

.19 

37 

SW-MA 

-.04 

10 

GE-AZ 

-.56 

10 

TC-NM 

-.35 

10 

GI-MA 

-.08 

5 

TFO 

-.32 

37 

GV-TX 

.17 

6 

TF-CL 

-.24 

15 

IIB-OK 

.14 

6 

TK-WA 

-.48 

5 

IIE-TX 

.43 

4 

TN-CL 

.29 

11 

IIK-WY 

-.17 

4 

TS-NI) 

-.24 

3 

HLrlD 

-.64 

49 

UBO 

.35 

35 

IIN-ME 

.29 

27 

VN-UT 

.22 

12 

HR- AZ 

-.37 

8 

VT-OR 

-.11 

11 

HY-MA 

-.62 

3 

WJI2YK 

-.'35 

8 

JE-LA 

.42 

3 

WI-NV 

-.10 

42 

JP-AT 

-.46 

8 

WMO 

-.01 

30 

JR-AZ 

.17 

10 

WM-AZ 

.20 

15 

KC-MO 

.15 

9 

WN-SD 

.16 

9 

KG-AZ 

.51 

4 

WO-AZ 

-.48 

9 

KM-CL 

.33 

4 

Ericsson  (1971)  has  obtained  the  same  result  with  magni¬ 
tudes  from  Canadian  stations  for  NTS  explosions.  Since 

/N 

our  C.’s  operate  somewhat  as  distance-effect  terms,  the 
interpretation  of  the  correlation  coefficient  is  tenuous; 
and  if  compressional-wave  station  terms  had  been  calcu¬ 
lated  relative  to  a  curve  or  to  line  segments  which  more 
closely  fit  the  actual  amplitude-distance  data,  the 
correlation  coefficient  might  well  have  been  higher. 
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STANDARD  DEVIATION  FOR  STATION  Ms 


figure  4a„  Histograms  of  standard  deviations  of  Ms  estimates 
for  41  NTS  events  using  JMD  and  conventional  methods. 


standard  deviations  for  78  explosions  are  calculated  from 
JMD  and  conventional  m^  estimates  with  the  same  data  base 
(all  of  the  events  of  Table  IV  were  used).  The  conven¬ 
tional  estimates  are  made  using  Evernden’s  (1967) 
distance-correction  formulas  for  the  Western  United  States 
and  Gutenberg  and  Richter’s  (1956)  teleseismic  distance- 
correction  terms.  Again  the  reduction  in  variance  is  due 
to  the  inclusion  of  station  terms  in  the  JMD  estimation. 
Ericsson  (1971)  shows  nearly  identical  results  using 
Canadian  magnitudes  for  NTS  explosions.  This  approximately 
100  per  cent  reduction  in  variance  for  Mg  and  m^  station 
estimates  should,  as  Ericsson  points  out,  correspondingly 
reduce  the  scatter  of  the  Mg  vs  m^  plot  using  network 
averages.  To  test  this,  we  took  36  of  the  43  explosions 
in  Table  I.  FISHER,  MADISON,  PAMPAS,  REDHOT,  and  STUTZ 
were  not  used  because,,  with  three  or  less  observations, 
their  M  was  not  well  estimated.  SEDAN  was  eliminated 
since  it  was  actually  a  cratering  event  and  appears 
anomalous  on  any  Mg  vs  m^  plot.  PAR  was  eliminated 
because  it  has  appeared  highly  anomalous  on  preliminary 
Mg  vs  m^  plots  (Lambert,  1971),  but  we  have  no  expla¬ 
nation  for  this  anomaly.  The  Rayleigh-wave  event  terms 
from  Table  I  and  corresponding  compress ional-wave  event 
terms  from  Table  IV  were  used  in  a  least-squares  pro¬ 
gram  which  fits  a  straight  line  to  data.  Note  that 
since  the  slope  is  nearly  one  for  the  fitted  lines, 
the  regression  of  m^  oh  Mg  could  not  differ  significantly 
from  that  of  Mg  on  m^  and  also  the  variance  of  m^  and  Mg 
about  the  fitted  line  would  be  nearly  equal.  Here  and 
in  all  subsequent  least-squares  fits  the  points  were  equally 
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weighted.  The  program  was  also  employed  to  fit  the 
“adjusted”  Mg  vs  m^  data  for  these  same  36  events,  as 
given  in  Lambert  (1971).  Using  the  truncated  data  sets 
of  this  report  (664  observations  for  Rayleigh  waves  and 
1356  observations  for  compressional  waves)  the  event 
magnitudes  were  computed  using  the  routine  formulas  just 
as  Lambert  did;  and  the  resulting  revised  Mg  and  for 
the  events  are  listed  in  Table  VI.  The  values  differ 
slightly,  or  not  at  all,  from  those  taken  from  Lambert; 
and  they  were  also  fitted  with  a  straight  line.  The 
fitting  results  are  listed  in  Table  VII  for  the  three 
cases;  and  the  actual  data  points  with  the  least-squares 
lines  are  shown  in  Figures  5,6,  and  7  for  Lambert’s 
values,  revised  Lambert’s  values,  and  JMD  values, 
respectively.  The  variances  about  the  lines  listed  in 
Table  VII  show  that,  some  difference  in  scatter  for  the 
three  plots  exists.  Little  decrease  in  variance  was 
obtained  by  using  magnitudes  estimated  from  the  truncated 
data,  which  comprised  about  93%  of  the  full  data  set  used 
by  Lambert.  Since  the  truncation  merely  removed  tails 
from  both  sides  of  the  distribution,  leaving  means 
essentially  unchanged  as  indicated  by  Table  VI,  little 
or  no  improvement  was  to  be  expected.  The  use  of  JMD 
results  versus  routine  magnitude  estimation  did  not 
decrease  the  variance  of  the  average  vs  m^  points, 
and  this  is  somewhat  surprising  in  view  of  the  approxi¬ 
mate  fourfold  decrease  of  the  variance  of  station  M„  and 

s 

m^  estimates  as  shown  by  Figures  4a  and  4b  w?*en  JMD  is 
employed. 

We  have  not  yet  shown  how  much  JMD  magnitude 
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TABLE  VI 


Comparions  of  Average  Magnitudes  for  36  NTS  Explosions 
When  Full  Data  Set  (Lambert)  and  Truncated  Data  Set 
(Revised)  are  used  with  Routine  Magnitude  Formulas. 


EVENT 

LAMBERT 

REVISED 

LAMBERT 

REVISED 

AUK 

3.97 

3.92 

4.50 

4.48 

BENHAM 

5.72 

5.79 

5.89 

6.28 

BILBY 

5.12 

5.12 

5.48 

5.45 

BOURBON 

4.49 

4.52 

4.78 

4.77 

BOXCAR 

5.86 

5.95 

6.14 

6.34 

BRONZE 

4.50 

4.51 

4,87 

4.86 

BUFF 

4.21 

4.22 

4.87 

4.86 

CHARCOAL 

4.00 

4.00 

4.82 

4.82 

CHARTREUSE 

4.33 

4.33 

5.04 

5.12 

COMMODORE 

5.24 

5.24 

5.40 

5,54 

CORDUROY 

4.77 

4.77 

5.42 

5.40 

CUP 

4.30 

4.29 

4.84 

4.83 

DILUTED  WATERS 

3.46 

3.50 

4.05 

4.05 

DUMONT 

4.69 

4.69 

5.27 

5.32 

DURYEA 

4.25 

4.21 

4.82 

4.P3 

FORE 

4.11 

4.09 

4.80- 

4.77 

GREELEY 

5.79 

5.83 

6.13 

6.17 

HARDHA.T 

*3.74 

3.72 

4.43 

4.34 

HAYMAKER 

3.98 

3.98 

4.34 

4.40 

HALFBEAK 

5.39 

5.36 

5.78 

5.93 

KLICKETAT 

4.04 

4.04 

4.57 

4.52 

KNICKERBOCKER 

4.78 

4.79 

5.22 

5.13 

MARSHMALLOW 

3.39 

3.27 

4.01 

4.02 

MERRIMEC 

3.45 

3.45 

3.73 

3.78 

MISSISSIPPI 

4.02 

3.99 

4.70 

4.75 

NASH 

4.14 

4.26 

4.93 

4.93 

PALANQUIN 

3.26 

3.26 

3.75 

3.76 

PILEDRJ.VER 

4.47 

4.47 

5.32 

5.33 

PINSTKIPE 

3.53 

3.49 

4.09 

4.11 

REX 

4.27 

4.26 

4.52 

4.55 

SCOTCH 

5,21 

5*13 

5.39 

5.58 

SCROLL 

3,06 

3.26 

3.89 

3.90 

TAN 

4.48 

4.49 

5.30 

5.28 

TURF 

4,06 

4,05 

4.55 

4.52 

WAGTAIL 

4.26 

4,26 

5.03 

5.04 

WISHBONE 

3.38 

3.42 

3.97 

4.01 

*This  value  differs  from  that  (3.88)  reported  by  Lambert  (1971), 
who  did  not  use  all  the  available  data.  This  value  would  have  resulted 
had  he  done  so. 
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Figure  5.  Ms  vs  mb  of  36  NTS  events  using  routine  magnitude 
estimation  --  from  Lambert  (1971). 
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Figure  6.  Ms  vs  mh  of  36  NTS  events  using  routine  magnitude 
estimation  with  the  truncated  data  sets  of  this  report. 


TABLE  VII 


Mg  vs  mb  Relations  Using  Various  Means  of 
Magnitude  Estimation  for  36  NTS  Events 


ESTIMATION 

_ PROCEDURE _ 

Routine  (Lambert's  original) 
Routine  (Lambert's  revised) 
JMD(A) 

JMD(A) 

JMD (A/T) 


M  vs  m.  RELATION 

S  D 

M„  =  1.08  m.  -  0.93 

S  D 

M_  =  1.03  m.  -  0.70 

S  D 

Ms  =  1.17  mb  -  0.60 

M_  =  1.41  m.  -  0.65 

s  b 

Mg  =  1.19  mb  -  0.57 


VARIANCE  OF  Ms  ABOUT 
THE  PREDICTED  RELATION 

.0395 

.0297 

.0366 

.0321 


0337 


estimates  differed  from  those  made  with  routine  formulas. 
To  do  this,  we  must  first  raise  the  36  event  terms  F  to 
a  level  equivalent  to  our  routine  and  mb  values.  We 
will  add  to  the  event  terms  a  term  R  given  by 

i  36  z\ 

'  Rs  ■  w  '  F)j 


for  surface-wave  magnitudes  and  by 
1  36  /s 

Rb  a  T5  X  Cmb  "  F)j 
J=1  J 

for  body-wave  magnitudes,  where  Me  and  m,  refer  to  the 

revised  values  given  in  Table  VI  and  F  refers  to  the 

JMD  event  terms  of  Tables  I  and  IV.  The  calculated 

value  of  Rg  was  4.34;  that  of  Rfe  was  4.39.  We  list  in 

Table  VIII  the  differences  in  routine  and  JMD  magnitude 

estimates  after  these  normalizing  R  terms  are  added  to 

/•> 

the  event  terms  F.  It  is  evident  that,  regardless  of 
the  fourfold  decrease  of  station  magnitude  variance 
with  JMD  due  to  station  terms,  the  JMD  average  event 
magnitudes  have  changed  very  little  from  conventional 
ones.  This  seems  to  account  for  the  failure  of  JMD  to 
provide  a  Mg  vs  mb  plot  with  less  scatter.  We  must 
conclude  at  this  point  that  the  observed  Mg  vs  mb  points 
in  Figures  5  and  6  are  real  and  relatively  unbiased  by 
the  particular  network  recording  each  event.  Since  JMD 
results  in  Figure  7  fail  to  reduce  the  scatter,  we  must 
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TABLE  VIII 


Difference  Between  Magnitudes  Computed  Routinely  and 
by  JMD  For  36  NTS  Events 


EVENT 

ROUTINE  - 
Ms-Fs 

JMD  ESI 

mb* 

AUK 

.03 

- .  06 

BENHAM 

-.15 

.16 

BILBY 

-.06 

-.02 

BOURBON 

,01 

.07 

BOXCAR 

-.10 

.23 

BRONZE 

.02 

-.01 

BUFF 

-.04 

-.05 

CHARCOAL 

.00 

.10 

CHARTREUSE 

-.05 

.00 

COMMODORE 

.04 

.02 

CORDUROY 

-.05 

-.05 

CUP 

.02 

-.03 

DILUTED  WATERS 

.12 

-.05 

DUMONT 

-.06 

-.01 

DURYEA 

-.03 

-.05 

FORE 

-.03 

-.10 

GREELEY 

.00 

.03 

HARDHAT 

.01 

-.03 

HAYMAKER 

.05 

-.14 

HALFBEAK 

-.09 

.03 

KLICKETAT 

- .  06 

-.13 

KNICKERBOCKER 

-.06 

.08 

MARSHMALLOW 

.02 

-.09 

MERRIMEC 

.07 

-.09 

MISSISSIPPI 

-.03 

- .  16 

NASH 

.05 

.02 

PALANQUIN 

.13 

.03 

PILEDRIVER 

-.10 

-.04 

PINSTRIPE 

-.03 

-.13 

REX 

.04 

.00 

SCOTCH 

.01 

.08 

SCROLL 

.15 

-.01 

TAN 

-.09 

-.08 

TURF 

-.01 

-.06 

WAGTAIL 

.02 

.03 

WISHBONE 

.12 

.08 

accept  such  scatter  and  attribute  it  to  differences  in 
source  parameters  rather  than  network  biases.  These 
parameters  are  source  depth,  medium  (including  density, 
rigidity,  bulk  modulus,  porosity,  water  content,  etc.), 
and  relative  tectonic  component;  but  it  is  beyond  the  • 
scope  of  this  report  to  relate  these  to  observed  magni¬ 
tudes  . 

Finally,  we  use  the  JMD  results  with  the  A  and  A/T 

measure  as  discussed  earlier  to  construct  Mg  vs  m^  plots 

in  Figures  8  and  9  respectively.  The  actual  F  terms  are 

given  in  Table  IX.  Straight-line  fits  were  made  to  these 

data,  and  the  results  are  listed  in  Table  VII  along  with 

the  other  M  vs  m,  lines.  Calculated  variances  indicate 
s  b 

that  the  scatter  of  these  last  two  plots  is  somewhat 
less  than  the  previous  one  where  A  was  used  as  the 
measure,  but  whether  this  is  significantly  less  is 
questionable* 
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Figure  8.  Ms  vs  mjj  of  36  NTS  events  using  JMD  results  with  A 
measure. 


Figure  9,  Ms  vs  mjj  of  36  NTS  events  using  JMD  results  with  A/T 
measure. 


TABLE  IX 


Event  Terms  From  JMD  Using  A  and  A/T  Measures 


AUK 

-  .46 

BENHAM 

1.70 

BILBY 

.86 

BOURBON 

.14 

BOXCAR 

1.72 

BRONZE 

.16 

BUFF 

-  .08 

CHARCOAL 

-  .29 

CHARTREUSE 

.03 

COMMODORE 

.82 

CORDUROY 

.48 

CUP 

-  .05 

DILUTED  WATERS 

-1.15 

DUMONT 

.41 

DURYEA 

-  .16 

FORE 

-  .18 

GREELEY 

1.50 

HARDHAT 

-  ,65 

HAYMAKER 

-  .41 

HALFBEAK 

1.14 

KLICKETAT 

-  .23 

KNICKERBOCKER 

.51 

MARSHMALLOW 

-1.18 

merrimec 

-1.09 

MISSISSIPPI 

-  .26 

NASH 

-  ,13 

PALANQUIN 

-1.37 

PILEDRIVER 

.28 

PINSTRIPE 

-  .87 

REX 

-  .17 

SCOTCH 

.80 

SCROLL 

-1,47 

TAN 

.25 

TURF 

-  .26 

WAGTAIL 

-  .06 

WISHBONE 

-1.36 

.44 
1.57 
.84 
.19 
1. 70 
.15 
••  .08 

-  .36 
.04 
.88 
.48 

-  .07 

-  .88 

.41 

-  .07 

-  .24 
1.40 

-  .61 

-  .41 

1.10 

-  .25 
.51 

-1.06 

-  .91 

-  .34 

-  .13 
-I. IS 

.21 

-  .80 

-  .10 

.85 

-1.15 

.23 

-  .30 

-  .12 

.93 


A 

.11 

1.49 
.96 
.33 

1.58 

.46 

.40 

.28 

.67 

1.00 

.90 

.42 

-  .31 
.85 
.41 
.35 

1.50 
.06 
.OS 

1.33 

.23 

.61 

-  .24 

-  .43 
.36 
.46 

-  .63 
.90 

-  .16 
.14 

1.08 

-  .43 
.75 
.10 
.58 

-  .31 


A/T 

.12 

1.64 

1.01 

.36 

1.66 

.46 

.47 

.30 

.73 

1.10 

1.02 

.43 

-  .33 
.91 
.  48. 
.40 

1.67 

.00 

.07 

1.45 

.23 

.68 

-  .29 

-  .51 
.40 
.52 

-  .61 
.97 

-  .13 

.20 

1.13 

-  .46 
.90 
.14 
.60 

-  .31 


ANALYSIS  OF  SCATTER  IN  Mg  vs  PLOTS 


We  desire  to  show  how  the  foregoing  results  fit 
into  an  amplified  statistical  discussion  of  M  and  m^. 
The  statistical  model  and  stated  results  in  this  section 
closely  parallel  the  work  of  Ericsson  (1971).  One 
important  distinction  is  that  he  used  known  yields  to 
estimate  the  variance  of  the  error  terms  in  the  model 
whereas  we  ignore  yield  information  and  obtain  nearly 
identical  results  from  the  seismic  measurements  them¬ 
selves.  Also  we  employ  a  significantly  larger  set  of 
explosions  (36)  in  the  analysis  of  variance  than 
Ericsson,  who  used  only  six.  We  will  compare  the  vari¬ 
ance.  of  single-station  Mg  vs  m^  plots  with  network 
Mg  vs  m^  plots  because  this  has  practical  importance 
and  was  not  discussed  by  Ericsson. 


First  we  take  each  station's  calculated  magnitude 
for  a  given  event  to  be  given  by  equation  (10) .  Then  we 
consider  Z..  to  be  composed  of  the  predicted  log-amplitude 

A  tj  « 

term  X..  plus  a  random  normal  error  term  E- .  peculiar  to 

a.  J  J 

that  station-event  pair  as  indicated  by  equation  (1).  Thus 


BR.  . 
ij 


H 


A 

ci 


(Ua) 


when  station  corrections  are  included;  or 


*  BRij  -  H  +  Btj 


(lib) 


when  they  are  not.  The  prime  henceforth  denotes  an  esti¬ 
mate  made  without  the  use  of  station  terms  C^.  E^.  is 
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"measurement  noise"  introduced  by  seismic  background  on 

the  recordings,  by  calibration  errors,  and  by  the 

analyst’s  subjective  measurement  of  amplitude  and  period 

/\ 

of  the  maximum  excursion.  is  "path  station  noise" 
which  is  due  to  varying  propagation  and  site  factors 
and  inability  of  a  single  straight  line  to  fit  the 
amplitude-distance  data  as  well  as  line  segments  or 
polynomials.  Applying  the  definition  of  to  (11a)  and 
(lib),  we  obtain  simply: 


F.  +  E.  . 
J  ij 


(12a) 


FL.  =  F.  +  E.  .  +  C.  (12b) 

J  11  i 


✓N 

The  Ej.  are  equivalent  to  those  in  equation  (1). 

Equations  (12a)  and  (12b)  relate  single-station 
magnitude  to  network  magnitudes.  Let  us  rewrite  them  as: 


AAA 


(13a) 


A  A  A  A 


(13b) 


We  now  consider  the  magnitude  differences  above  as  ran¬ 
dom  variables  with  variances  given  by: 
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9  A  A  O  A 

9  A  A  9  A  9  A 

*i(Fij  -  V  -  ci<V  +  ai<V 


2 

where  the  subscript  i  on  a  indicates  a  summation  over 
stations.  From  the  final  iteration  of  the  JMD  program, 

A 

we  obtain  the  variance  of  all  the  E..’s,  and  the  vari- 

A  1J 

ance  of  the  C^'s  is  easily  obtained  from  the  values  of 

these  terms  as  given  in  Tables  II  and  V.  The  computed 
2  A 

values  of  o..(E..)  are  .020  for  surface  waves  and  .036 
i]  2  A 

for  body  waves,  and  the  values  of  o.(C^)  are  ^ 046  for 

surface  waves  and  .150  for  body  waves.  (For. the  body 
2  A 

waves,  o^(C^)  should  be .significantly  reduced  by  using 
line  segments  or  a  polynomial  to  fit  the  amplitude- 
distance  data.)  These  values  are  reflected  in  Figures  la 
and  lb,  2a  and  2b,  4a,  and  4b  which  compare  scatter 
before  and  after  station  corrections  are  applied  in 
magnitude  estimation. 


The  can  also  be  relate.d  to  the  scatter 

observed  in  Mg  vs  m^  plots  for  single  stations  compared 
to  those  for  networks  as  in  Figures  10a,  10b,  and  10c. 
In  each  figure  the  identical  event  set  is  used  for  both 
single-station  and  network  estimates  of  magnitude;  the 
events  and  estimates  are  listed  in  Table  X.  The  scatter 
in  the  single-station  plots  is  greater,  and  we  can 
express  this  quantitatively  by  forming  equation  (12a) 
for  both  body  and  surface  waves  and  subtracting  to 
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Figure  10b.  Ms  vs  mj,  at  NP-NT  compared  to  network  estimates 
for  same  events. 


TABLE  X 


Event  Terms  C.  ulated  from  Three  Single 
Stations  and  From  Entire  Network 


RK 

Fs 

-ON 

Fb 

HN 

FS 

-ME 

Fb 

NP' 

Fs 

-NT 

Fb 

NETWORK* 
Fs  Fb 

BENHAM 

1. 38+ 

1.49 

1.60 

1.39 

1.60 

1.73 

BILBY 

.80 

1.07 

- 

- 

.81 

.92 

.84 

1.08 

BOURBON 

•» 

•» 

-  .05 

.61 

.24 

.25 

.17 

.31 

BOXCAR 

- 

- 

1.85 

1.43 

2.04 

2.11 

1.71 

1.72 

BRONZE 

.24 

.24 

.10 

.85 

.18 

.42 

.15 

.48 

BUFF 

.11 

.34 

-  .31 

.66 

-  .21 

.40 

-  .08 

.52 

CHARCOAL 

.31 

.14 

-  .SI 

.61 

- 

- 

-  .34 

.33 

CHARTREUSE 

- 

- 

.17 

.55 

-  .07 

.78 

.04 

.73 

COMMODORE 

.91 

1.09 

.83 

.82 

.97 

1.14 

.86 

1.13 

CORDUROY 

.51 

1.01 

- 

- 

.62 

1.14 

.48 

1.06 

CUP 

.11 

.43 

-  .21 

.45 

-  .18 

.40 

-  .07 

.47 

DUMONT 

.68 

.93 

.44 

1.21 

.37 

.70 

.41 

.94 

DURYEA 

.14 

.89 

-  .26 

.13 

-  .12 

.87 

-  .10 

.49 

GREELEY 

- 

- 

1.59 

1.57 

- 

- 

1..49 

1.75 

HALFBEAK 

- 

- 

1.31 

1.14 

- 

- 

1.11 

1.51 

KNICKERBOCKER 

.44 

.95 

.69 

.81 

.68 

.87 

.51 

.66 

MISSISSIPPI 

- 

- 

-  .21 

.47 

- 

- 

-  .32 

.52 

NASH 

- 

- 

- 

- 

-  .22 

.32 

-  .13 

.52 

PILEDRIVER  - 

.08 

.65 

.54 

.96 

.07 

.69 

.23 

.98 

SCOTCH 

.98 

1.53 

.93 

1.11 

.80 

1.39 

.83 

1.11 

TAN 

.43 

1.05 

1.09 

1.18 

.00 

.67 

.24 

.97 

TURF 

- 

- 

-  .49 

-  .04 

- 

- 

-  .28 

.19 

+A11  F  estimates  include  terms  (s). 

*Network  estimates  are  same  as  given  in  Tables  I  and  IV. 


-Lacks  compressional-wave  or  Rayleigh-wave  recording,  or  both. 


obtain: 


CF3  -  Pj) 


AC 

+  E.  . 
13 


Considering  the  above  magnitude  difference  as  a  random 
variable  again,  we  can  express  the  variance  as: 


_2 ,£s  £b  s  _  rr2r~s  ^bs  .  2,£s  ,  .  n2*b  , 

*  Fij)  -  VF:  '  Fj5  +  °j  CEij5  +  W 


2 

where  the  subscript  j  on  a  indicates  a  summation  over 
events.  Using  the  computed  variances  of  the  magnitude 
differences  as  given  in  Figures  10a,  10b,  and  10c,  we 
obtain  Oj(E?j)  +  o?(E^)  equal  to  .039,  .033  and  .070 
for  RK-ON,  NP-NT  and  HN-ME  respectively.  (The  computed 
variances  on  these  figures  are  for  the  vertical  devia¬ 
tion  of  the  points  from  the  least-squares  fitted  lines. 
(Using  the  vertical  distance  actually  gives  the  variance 
of  Mg  -  a  m^  where  a  is  the  slope  of  the  fitted  line. 
None  of  the  analysis  is  seriously  affected  by  values 
of  a  which  are  slightly  different  from  1.0.)  We  know 
from  the  JMD  results  that  a?.  (E?.)  +  cr?.  (E^.)  is  .056 
for  the  entire  group  of  data  covering  more  than  a 
hundred  stations.  Thus  RK-ON  and  NP-NT  exhibit  less 
measurement  error  than  at  a  typical  station  while  HN-ME 
exhibits  more. 
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From  equation  (12a)  we  have  shown  that  scatter  for 
single-station  M0  vs  mw  plots  will  exceed  that  for  a 
network  byaj(E^.j)  +  o?(E_)  f0r  that  particular  station 
when  station  corrections  are  applied.  When  they  are  not, 
we  must  proceed  from  equation  (12b)  and  let 


A 


Li 

l 5  C. 

i=l  1 


(14) 


where  Lj  is  the  number  of  stations  recording  the  j’th 
event.  Substituting  (14)  in  (12b),  forming  (12b)  for 
body  waves  and  surface  waves  again,  and  subtracting,  we 
obtain: 


i 

(Ft . 


Ft  .) 
ir 


As» 


l! 


CFJ 


Lj 


1_  T3?S  ib'  1 

7J  .Lci  Fj  +  7^  ,LL i> 

Lj  1-1  Lj  1-1 


As  Ab  Ac  Ah 

*  EIj  -  Eij +  ci :  ci 
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We  can  write  the  scatter  of  this  magnitude  difference 
as : 


*s» 
o  .  (F .  . 

r  13 


PbS  -  rt2r5s»  nb’ 

V  '  VFj  '  Fj 

Lb 

*  h  r  £ i) 

L?  i=i  1 


L*  i=1 


.2,5s 


+  + 


.2,5s 


«jfcp  ♦ 


(15) 


A 

Since  C.  is  constant  from  event  to  event,  its  variance 

X  A 

is  zero  here.  By  equation  (14),  Fj  is  defined  as  the 
sum  of  two  random  variables,  and  its  scatter  must  be 

A 

greater  than  that  of  Fj  so  that 
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Applying  the  above  principles  to  (15) ,  we  have 


2/cl  Ahl  ■>  A  | 


AK  »  O  A 

-  FT  )  +  a r(E?0 
3  3K  13 


CEij) 


3  L!  i=1 


Using  the  fact  that  the  variance  of  the  mean  of  L. 
realizations  of  a  random  variable  is  equal  to  1/L^ 
times  the  variance  of  the  variable  itself,  this  reduces 
to: 


a?(F?!  -  F1?!)  =  a?(F?f  -  Fb,)+  a2(Es  ) 
r  13  i3J  V 3  3  jU ir 


(16) 


L; 


L. 

3 


This  relation  shows  that,  when  station  corrections  are 
not  applied  and  the  network  has  a  large  number  of  sta¬ 
tions,  the  scatter  of  a  single-station  Mg  vs  mb  plot 
will  exceed  that  of  a  network  by  the  combined  variance 
of  the  measurement  error  as  in  the  case  when  station 
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corrections  are  applied.  However,  if  Lj  or  Lj  is  small 
and  station  corrections  are  relatively  large  compared  to 
measurement  errors  the  single-station  plot  may  actually 
have  less  scatter  than  the  network  plot. 


We  desire  to  add  still  another  noise  term  to  our 

/N 

analysis  by  letting  F.  be  a  linear  function  of  the  log- 

J  A 

arithm  of  yield  plus  a  random  variable  Sj  which  is  the 
aberration  from  the  expected  magnitude  for  a  given 
yield  due  to  detonation  depth,  source  medium,  and 
possibly  tectonic  strain-release  influences.  For  a 

A 

real  event,  there  may  be  some  variation  of  the  among 
stations  due.  to  a  dependency  of  this  factor  on  azimuth 
and  take-off  angle  of  the  ray  path  at  the  source,  but 
we  will  disregard  such  variation  in  this  discussion. 

'  '  A 

In  fact,  if  such  variation  in  S.  is  consistent  amon' 

J  A  A 

events,  it  will  be  absorbed  by  the  terms.  The 
are  estimatible  only  if  the  yields  are  known;  however, 

A 

we  can  estimate  the  variance;  of  the  S^'s  without  yield 
information.  The  linear  relation  between  magnitude  and 
yield  has  a  substantial  theoretical  ar.<\  empirical  basis 
up  to  approximately  100  kt  or  m^  ~  6  (von  Seggern  and 
Lambert,  1970;  Lieberroann  and  Pomeroy,  1969;  Evernden, 
1969;  Evernden  and  Filson,  1971;  Ericsson,  1971).  From 
these  studies,  the  coefficient  on  the  logarithm  of  yield 
ranges  from  approximately  0.8  to  1.2.  Without  any  sig¬ 
nificant  effect  on  our  results  we  can  assume  it  is  1.0 
for  both  body  and  surface  waves.  We  make  these  defini¬ 
tions  then: 


A 


log  Yj  + 


K  +  +  rr 

J  J  * 


Li 

J  A 

l  E.  . 
i=i  ^ 


(17a) 
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0 


FI 

3 


Li 

3* 


L, 

JA 


-  log  Y.  +  K  +  S.  +  l  E.  .  +  r—  l  C. 
3  3  Lj  13  L4  4-i  1 


3  i-1 


These  definitions  include  the  mean  of  the  measurement 

A 

error  terms  since  the  JMD  estimate  of  Fj  will  neces¬ 
sarily  be  biased  by  this  value^  and  we  do  not  want  to 
include  it  in  the  source  term  Sj.  If  we  form  (17a)  and 
(17b)  for  both  body  and  surface  waves  and  subtract,,  we 
obtain: 


*  kS 


K  = 


As  Ah 

sj  •  si  + 


Lb 

L.  sV. 
If  i=l  ^ 


i_  j3e*. 

LS  i=l  x3 
3 


(18a) 


As  *  Ab » 
F.  -  F . 


+  KS  -  Kb  =  s! 


Lb 


L® 

Sb  +  ^E? . 

j  Lj  iii  V 


L? 


Lb 


ij-  .  +  i-  -  ir~  ^Cb 

ib  ij  s  i  7b  .  i 

x  *”*  X  L  •  1  ~  1  L  *  1 w  X 


(18b) 
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The  scatter  of  these  magnitude  differences  is: 


2,£s  Cbs  _  _2,£s-v  .  _2,$bv  1  2,~s  1 

VFj  "  V  "  W  a3CS3:>  +  7T 


u  1,  2r£b  s 

Lj 


(19a) 


rt2,Cs»  CbfN  _  „2fZs,  .  „2,Zb,  .  1  2,£s  > 
aj(Fj  -  F.  )  -  0  j  (Sj)  +  cr.CSj)  +  -  ^(E..) 

j 


(19b) 


+  7  ai(Eii}  +  ? 
j  1 


2  rpS  -v 

°i  (Ci' 


.  1  rt2,pb. 
rh  ai^ci) 
Lj 


where  again,  as  in  arriving  at  (16) ,  we  used  the  fact 
that  the  variance  of  the  mean  is  1/Lj  times  the  variance 
of  the  variable  itself. 

We  can  now  estimate  the  variance  of  the  source  terms 
by  applying  equation  (19a)  to  the  JMD  results  shown  in 
Figure  7.  There  the  variance  of  the  points  about  the 


-32- 


fitted  line,  measured  vertically,  is  approximately  .037. 
Variances  on  other  Mg  vs  m^  plots  for  these  same  events 
in  Figures  5,  6,  S  and  9  are  within  20  percent  of  this 

c  K 

value.  On  the  average,  L?  =  17  and  L;:  =  20  for  the  36 

J  O  A.  1 

events  of  Figure  7;  and  using  ot(E. .)  =  .020  and 

2  Ab  i  ij 

°i^Eij^  =  *936  *n  ecluation  (19a)  >  we  estimate  the 

combined  source,  variance  to  be  .034.  Since  we  have  no 

means  of  separating  the  two,  we  simply  assume 

a? (Sp  =  0?  (Sp  =  .017.  With  yield  information  Ericsson 
j  j  J  J  o  As 

was  able  to  make  the  separation  and  obtained  of( Sp  =  .012 

2  ,Ab  J  J 

and  -  .017  for  his  data.  It  is  noteworthy  that 

the  six  explosions  which  Ericsson  used  were  all  detonated 
within  a  small  test  area  on  Pahute  Mesa  in  either  a 
rhyolite  or  tuff  medium  and  that  the  36  explosions  con¬ 
sidered  in  this  report  were  detonated  over  a  much  greater 

area  and  in  a  variety  of  media.  On  that  basis,  we  expect 
2 

cK(Sj)  to  be  much  larger  for  our  group,  hut  it  is  in 
fact  only  slightly  so. 

We  now  can  explain  the  failure  of  station  correc¬ 
tions  to  improve  the  network  Mg  vs  m^  scatter  as  shown 
by  Figures  6  and  7.  Equation  (19b)  applies  to  Figure  6, 
while  equation  (19a)  applies  to  Figure  7,  Using 

O  A _  o  AL 

07  (Cp  =  .046  and  07  (Cp  a  .150  as  previously  determined, 
we  see  that  since  these  are  divided  by  the  number  of 
stations,  their  combined  weight  in  equation  (19b)  is 
about  one-third  to  one-fourth  the  combined  weight  of 

9  Ac  9 

o.(Sp  and  0-(S.),  with  the  measurement  variances  also 
j  j  J  J 

relatively  small.  It  is  evident,  though,  that  their 
weight  will  increase  as  the  number  of  stations  decreases, 
so  that  station  corrections  are  imperative  when  small 
changing  networks  or  random  isolated  stations  record  events. 
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CONCLUSIONS 


The  JMD  method  of  magnitude  estimation  is  ,a  simple 
and  straightforward  procedure  which  performs  as  well  as 
routine  magnitude  estimation.  This  is  shown  by  the  fact 
th„c  Mg  vs  m^  plots  show  about  the  same  scatter  when  JMD 
results  for  36  NTS  explosions  were  compared  against  rou¬ 
tine  ones.  Because  of  its  manner  of  incorporating  distance 
effects  into  station  corrections,  th®  method  can  only  be 
applied  uniformly  to  a  limited  source  region,  perhaps 
only  a  hundred  kilometers  wide. 

Both  Rayleigh-wave  and  compressional-wave  ampli~ 
tudes  from  NTS  explosions  were  shown  to  be  log-normally 
distributed  over  a  large  station  network.  This  fact 
considerably  simplifies  the  statistical  model  for 
magnitudes  and  justifies  the  traditional  procedure  of 
averaging  station  magnitudes  to  obtain  an  event  magnitude 
estimate. 


In  our  analysis  of  variance  for  Mg  and  m^, 


var (surface-wave  measurement  noise)  =  .020 
var (body-wave  measurement  noise)  =  .036 
var (surface-wave  path  station  noise)=  .046 
var (body-wave  path  station  noise)  =  .150 
var (surface-wave  source  noise)  =  .017 
var (body-wave  source  noise)  =•  .017 


we  found: 

=  var(EijS) 
=  var(E„b) 
=  varfC^) 

=  varfC^k) 

=  var (S.s) 

A  U 

=  vav(SjD) 
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.One  can  estimate  these  variances  without  the  knowledge 
o£  the  explosive  yields,  provided  there  is  a  sufficient 
number  of  events  and  stations  to  apply  statistical 
procedures  reliably.  We  emphasize  that  these  results 
apply  only  to  the  data  base  used  here  --  LRSM  stations 
and  NTS  explosions  --  and  may  change  significantly  with 
other  data  bases. 

Of  the  three  types  of  noise,  the  station-path  term 

A 

is  predominant.  Its  biasing  of  magnitude  can  be 
eliminated  entirely  by  repeated  use  of  a  single  station 
or  of  an  identical  group  of  stations  to  estimate  event 
magnitudes;  in  this  case  all  magnitudes  will  be  biased 
by  an  equal  amount  and  thus  not  contribute  to  scatter 
in  Ms  vs  m^.  Its  biasing  effect  can  also  be  suppressed 
by  use  of  well-determined  station  corrections.  If  the 

/N 

term  can  be  kept  small  by  accurate  recording  system 

calibration  and  careful  systematic  measurements,  even  a 

single  station  could  give  Mg  vs  m^  plots  with  a  small 

scatter.  This  can  be  observed  for  Berkeley  recordings 

(McEvilly,  1971)  and  Lamont  recordings  (Savino,  1971), 

When  compelled  to  use  varying  stations  or  networks  for 

magnitude  estimation,  as  many  stations  as  possible 

✓\ 

should  be  used  because  the  effect  of  the  C-  term  is 

1 

reduced  by  increasing  the  number  of  stations.  This  holds 

A 

for  the  E^  term  also.  If  a  large  number  of  stations  are 
used,  say  ten  or  more,  then  the  application  of  station 
corrections  clearly  becomes  superfluous  in  reducing  Mg 
vs  m^  scatter. 

/V 

The  source-term  variance  Sj  cannot  be  suppressed  and 
is  the  primary  remaining  cause  of  Mg  vs  m^  scatter.  This 
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alone  gives  a  standard  deviation  of  0.13  units  normal  to 
the  Mg  vs  m^  straight-line  fit  for  the  NTS  data  (assum¬ 
ing  a  slope  of  unity).  We  feel  that  this  is  less  than 
the  scatter  that  could  result  if  a  conscious  continuing 
attempt  were  made  to  select  the  detonation  media  and 
depth  in  such  a  way  as  to  vary  the  source  coupling  in 
short  and  long  periods.  An  examination  of  source  medium 
and  depth  effects  on  Mg  and  m^  coupling  seems  essential 
if  we  are  to  define  the  expected  limits  of  possible 
scatter  in  explosion  Mg  vs  m^  plots  for  a  selected  test 
region. 
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